Specific features of the effect of time dependent field on subdiffusing particles. The 

stochastic Liouville equation approach. 



OO 
O 

o 

(N 

ojq: 
< 

(N 
(N 



i 

c 

o 
o 



o 

cn 

o 

OO 
O 



X 



A. I. Shushin 

Institute of Chemical Physics, Russian Academy of Sciences, 
117977, GSP-1, Kosygm str. 4, Moscow, Russia 

We analyze the effect of time dependent external field on non-Markovian migration described 
by the continuous time random walk (CTRW) approach. The rigorous method of treating the 
problem is proposed which is based on the Markovian representations of the CTRW approach and 
field modulation. The method is applied to the case of subdiffusive migration in which the exact 
formulas for the first and second moments of spatial distribution are derived. For oscillating external 
field they predict unusual dependence of the first moment on oscillation phase and anomalous field 
dependent contribution to the dispersion. Similar formulas are also derived fluctuating field. 

PACS numbers: 05.40.Fb, 05.40.Jc, 02.50.-r, 76.20.+q 



I. INTRODUCTION. 

Brownian motion in external time-dependent field is 
the important stage of many physical and chemical pro- 
cesses which often strongly affect their kinetics 0, 
Last years close attention is given to the anomalous 
(subddifusiye) jump-like motion typical for disordered 
systems [!, Q and, in particular, to the effect of time- 
dependent field on this type of migration dH]. Usually 
the motion anomaly is assumed to be a manifestation 
of the long memory in the kinetics of jumps. In such a 
case the serious difficulty in theoretical treatment of time- 
dependent field effects occurs because of subtle interplay 
of field and anomalous memory effects which should be 
properly described. 

Subdiffusive processes in time-independent potential 
V(x) are traditionally described by the fractional Smolu- 
chowski equation (FSE) for the probability distribution 
function (PDF) p(x, t) [| 



P=~(Pt ° C <*P> 



(1) 



where gD t a is the Riemann-Liouville fractional deriva- 
tive defined by 



d 



and 



T(a)dtJ '(t-ti) 1 -" 



-D a V x [V x - F(x)} 



(2) 



(3) 



is the Smoluchowski operator, in which D a is a subdiffu- 
sion constant, V x = d/dx, and = — \7 x V(x)/(kBT) 
is a force. The FSE ((TJ can be derived within the contin- 
uous time random walk (CTRW) approach 0] assuming 
the long time tailed behavior of the waiting time distri- 
bution W(t) for CTRW-jumps: W(t) ~ l/t 1+a (a < 1). 

In the case of time-dependent field F(x, t) [i.e. time- 
dependent C a (i)), however, no analogs of the FSE are 
rigorously derived as yet. The main difficulty is in correct 
treatment of the effect of £ a (i)-evolution during the time 



of waiting for jumps. Only approximate variants of these 
FSEs have been proposed so far [f| [|| . 

In this work within the CTRW approach we derive the 
exact FSE (describing the influence of time-dependent 
field) with the use of recently proposed Markovian repre- 
sentation of the CTRW and the non-Markovian stochas- 
tic Liouville equation (SLE) [7]. The solutions of this 
FSE for different time dependences of the force F(t), for 
simplicity assumed to be independent of x, are proposed 
and discussed in detail. 



II. MARKOVIAN SLE. 

Here we will present the method of treating the effect 
of time-dependent field F(x, t) on CTRW-like processes 
by reduction of the problem to solving the SLE with time- 
independent operators. 

To clarify the method we first consider the Markovian 
(normal diffusion) case: a = 1, in which the evolution of 
the system is described by the Smoluchowski equation 



P : 



-Ci(t)p = D 1 V x [V x p - F(x,t)p]. 



(4) 



The method is based on the representation of the time 
dependence of F(x,t) in terms of the dependence on 
some Markovian (in general, stochastic) variable z(t): 
F(x,t) = F(x,z(t)), whose evolution is described the 
PDF a(z,t) satisfying the Markovian equation 



-La, with a(z, 0) = <xs(z), 



(5) 



in which L is the linear operator in {z}-space and 
J dztji(z) = 1. For brevity, formulas are written assum- 
ing that {z}-space is one dimensional, though they are, 
evidently, valid for any dimensionality of {z}-space. The 
corresponding examples will be discussed below. In ad- 
dition, in what follows we will restrict ourselves to the 
simple model of x-independent force: 



F(x,t) = F(x,z{t)) = F z(t). 



(6) 
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In this representation the kinetics of the process de- 
scribed by eq. ^ is determined by the average evolution 
operator which in the space {x ® z} is given by formula 

U(x,z;Xi, Zi \t) = (x,z\(Te-^ dT ^)\x t , Zl ), (7) 

where the average (denoted as (. . . )) is taken over trajec- 
tories of the stochastic Markovian process in {z}-space 
with fixed initial (zi) and final (z) coordinates. In par- 
ticular, the PDF of interest, p F (x,Xi\t) averaged over 
^(^-fluctuations, can be calculated as 



p F (x,Xi\t)= J dzj dziU(x,z;xi,Zi\i)<Ji(zi). (8) 

The important point of the proposed representation 
consists in the fact that for Markovian processes in {x (g) 
z}-space the operator U is satisfies the Markovian SLE 
with time-independent operators Q 



U = -(C 1 +L)U 



(9) 



with U(x, z; Xi, 2j|0) = 5(x — Xi)5(z — Zj). 

Thus we have reduced the problem to solving the SLE 
([9]) with time-independent operators, though at the cost 
of the extension of the space of the process, which de- 
scribes the evolution of the system. 

Noteworthy is that the representation ©-(JSJ) is valid 
not only for stochastic functions z(t) but also for dynam- 
ical ones, which are known to be Markovian as well. For 
example, in the model of harmonically oscillating force: 



z(t) = z c (t) = zq sm(ujt + ip), 



(10) 



the dependence z c (t) can be considered as a coordinate 
part of the trajectory of dynamical motion (in the har- 
monic potential), described by the operator 



L = -(uV z + lu 2 zV Vz ), 



(11) 



in the phase space {z} = (z,v z ) (v z = z is the velocity) 
with <7j(z) = 5(z — ZQsincp)5(v z — Zocostp). Evidently, the 
case of z(t) represented as a linear combination of, say, 
N oscillating functions Zj(t) — z . sm(ujjt + ipj): z(t) — 

Zj(t), can be modeled by coupling to N harmonic 
coordinates zjv = {z\, Z2, ■ ■ ■ , zn). 



III. NON-MARKOVIAN CTRW. 

A. Markovian representation. 

The main goal of this work is the analysis of the effect 
of time-dependent field on CTRW-type (subdiffusive) mi- 
gration. 

In the CTRW approach the stochastic motion in {x}- 
space is treated as a set of jumps with jump statistics 
described by the waiting time distribution W(t) [!, 0|. 
For time-independent driving force the non-Markovian 



equation for the PDF p(x, t) in is conventionally derived 
by summing up the contributions of all sets of jumps. 
In terms of the Laplace transform R(e) = f£°dt p(t)e~ et , 
this equation is written as [1, H[ 

eR(e) = Pi - M(e)£ a R(e). (12) 
In this equation pi(x) is the initial PDF and 



M(e) = [l-W(e)]/[eW(e)}, 



(13) 



where W(e) = J °° dtW(t)e~ et . Note that in the case of 
subdiffusion, when M(e) = e 1_Q , eq. (fT!?)) reduces to the 
Laplace transform variant of the FSE (TT|) . 

CTRW-type processes can conveniently be analyzed 
within the Markovian representation Q in which these 
processes are assumed to result from jump-like Ci{t)- 
fluctuations determined by the dependence of C\{y{t)) 
on some Markovian stochastic variable y{t) whose PDF 
i](y,t) satisfies equation 



f] = -Ai], with r)(y\Q) = rji(y). 



(14) 



In this equation A is a linear operator in {y}-space and 
r]i(y) is the initial condition [Jdyi]i(y) = 1]. 

The dependence C\{y) is taken in the form C\(y(t)) = 
5(yo—y(t))Cx, where yo is the coordinate at which the 
system undergoes the jump described by C\. Similar 
to the above-considered model of z(£)-modulation of Ci, 
in the case of j/(t)-modulation the evolution of the sys- 
tem is described by the PDF p(x,y\t) in the combined 
space {x <g) y}. This PDF satisfies the SLE of type 
of eq. ©, which as applied to the Laplace transform 
P(e) = f^°dtp(t)e~ et is given by 

eP = 5(x-Xi)5(y-Vi) - [A + 5{y-y )ti]P. (15) 
Of special interest is the PDF averaged over y(t)-process 

p y (x,Xi\t) = j dy JdyiPix^^Xi^tltfrniyt). (16) 

The y(t)-controlled (or modulated) process in {x}- 
space proves to be of CTRW type. Thus obtained CTRW 
depends on the initial condition pi (y) and the form of the 
operator A. In what follows we will consider the non- 
stationary variant realized for Pi(y) = S(y — yo) 0- l n 
this variant the average PDF ~p y (x, Xi\t) is known to sat- 
isfy the CTRW-like equation usually written as applied 
to R(e) = R v (e) = f™dtp- v (t)e-* 0: 

eR{e) = Pi - M(e)C a R(e), (17) 

in which 

M(e) = M(e) = (D,/ D a )e(y \(e + A)" 1 ^). (18) 

The behavior of M(e) is completely determined by the 
specific features of the controlling process (fT4]) . Various 
models of this process are discussed in ref. [7[ • 
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IV. CTRW-BASED SLE. 

The Markovian representation is very suitable for 
treating the effect of time-dependent field on CTRW-like 
processes. 

The corresponding equation is straightforwardly de- 
rived by taking into account that, in accordance with 
the Markovian representation, this equation describes 
Markovian process in {x ® y}-space affected by the the 
driving force which can be modeled by interaction with 
the Markovian z(t) variable. This means that the equa- 
tion sought is actually the Markovian SLE Q for the 
PDF q(r; r,;|i) in the extended space {r} = {x ® y <g> z}- 
space. For the Laplace transform Q(e) = J Q dtq(t)e~ et 
this equation is written as: 

eQ = Srr, - (A + Syy.tx + L)Q (19) 

This equation is seen to differ from cq. (fTS"]) for 
P{x, y\e) only in the replacement e by fl = e + L and the 
evident change of (S-function describing the initial condi- 
tion. Naturally, for the PDF averaged over y(t)-process: 

/■oo 

R{e)=R y {e)= / dtp y (t)e~ et , (20) 
Jo 

one gets CTRW-like equation (sometimes called the non- 
Markovian SLE 0]) similar to eq. (17]) 

[Cl + £ a M(tl)]R = pm with (l = e + L. (21) 

Notice that the order of operators C a and M(fi) is im- 
portant since they do not commute with each other. 

To qualitatively interpret eq. (f2Tj) within the CTRW 
approach it is worth noting that, according to the 
SLE representation (|9|) , the time-dependent-field affected 
CTRW can be considered as sequence of jumps governed 
by the z-dcpcndcnt operator C\{z). The CTRW pro- 
cess is accompanied by the simultaneous evolution of the 
parameter z(t), determined by the operator e~ Lt . This 
operator will enter in formulas describing CTRW evolu- 
tion in the form of the product W(t)e~ Lt which for the 
Laplace transform Wie) just corresponds to the replace- 
ment e by fl = e + L in W(e) thus resulting in eq. (|21|1 . 

Exact eq. (|2"Tj) essentially differs from that proposed 
earlier 0, Q to treat the effect of time-dependent field, 
i.e. the results of these works are, in general, incorrect, 
except, may be, some special cases (see below). 

A. General results. 

In this brief communication we will restrict ourselves to 
analyzing the time evolution of the moments of the PDF 
p yz (x\t) (averaged over y(t)- and z(i)-processes) m n (t) — 
J dx x n p yz (t) . For this analysis we need to specify the 
initial PDF pt(x). For simplicity we will assume pi{x) = 
5{x). 



Instead of moments m n (i), it is more convenient to 
analyze their Laplace transforms m n (e) which can be ex- 
pressed in terms of moment operators in the {z}-space 

M n = Sdxx n R~ y {x,(l) : (22) 

m„(e) = (M n ) z = Jdzfdz(z\M n \z)(7i(z). (23) 

As is seen from eq. (|2"Tj) the operators M n (e) satisfy 
simple equations 

QM n = n/^M(0)M„_i + n(n - l)D a M n - 2 (24) 

for n > 2 with M-\ = and Mo = O -1 , in which 
/ = D a FoZQ. Solution of these equations and substi- 
tution into eq. (|23l) yields for Laplace transforms of time 
derivatives of the moments 

riii(e) = /<z*(ft)>„ ^ 2 (e)=A (e)+A 2 (e), ( 25 ) 
with 

Ao(e) - 2D a {$(£l)) z , l 2 (e) = / 2 ((z$(0)) 2 },. (26) 

Here $(f2) = Cl^M (fi). 

In what follows we will concentrate on the analysis of 
force dependent terms. The force independent contribu- 
tion fi (e) was discussed in detail earlier [J. 

After the inverse Laplace transformation of expressions 
([251) and (261) one gets fi G {t) = 2D a <j)(t), 

mi(t) = fjjdzdzz zg(z,Zi\t)ai(zi), (27) 
P2{t) = f 2 JJJdzdzdzi zz 

x f*dig(z,z\t-F)g(z,z i \t)ai(zi). (28) 

In these formulas 

g(z,z\t)=<P(t)(z\e- Lt \z), (29) 

where (z\e~ Lt \z) is controlled by the model of z(t)- 
modulation, while <j>(t) = (2ni)^ 1 f^^de e~ 1 M(e) is de- 
termined by the CTRW model considered and for the 
subdiffusion model \M(e) — e 1_a , (a < 1)] 

<t>{t) =T-\a)t a - 1 . (30) 

B. Applications 

a. Harmonically oscillating force. In the model of 
harmonically oscillating force (flTJ|) in which L [eq. (|lip] 
describes dynamical motion in the harmonic potential, 
one gets (z|e _i *|z) = 5(z — z c (z\t)), where z c (z|t) = 
(z c (z\t),v Za (z\t)) is the trajectory of dynamical motion 
with zjt = 0) = z in the phase space {z}. 

Substituting this formula into eqs. (|27p - ([2"9"|) one ob- 
tains no(t) = 2D Q F- 1 (1 + a)t a , 

mi(t) = ff*dTZ c (r)4,(T), (31) 

p 2 (t) = f 2 f*diz c (t)f*dTZ c (T)<t>(i-T)<j>(T), (32) 
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where z c (t) is given by eq. (fTO|) . 

For brevity, we will restrict ourselves to the discussion 
of the long time behavior of the moments only: 

mi(t) (f/Lu a )sm((p + Tra/2) (33) 

m (t) ^ 72 (a)/ 2 (^) Q (34) 

where / = D a F zo and 72(a) = cos(7ra/2)/[2r(l + a)]. 

These formulas predict some peculiarities of the subd- 
iffusion response to oscillating force. First, mi(t) appears 
to be nonzero with asymptotic value (at t — > 00) indepen- 
dent of time and harmonically oscillating as a function of 
the initial phase ip of force oscillations. Second, the force 
dependent part of pi-zif) is anomalously large increasing 
in time so that [i 2 (t)/no (i) = cos(na/2)[f 2 / (4D a ui a )] 
independent of time. Third, in the case of conventional 
diffusion, i.e. at a — ► 1, ^i 2 (t)/t a — > 0, as should be. 

The exact formulas (|3"3"|) and (j3~4")l significantly differ in 
their analytical form from those derived earlier with not 
quite correct kinetic equation [5j], though, surprisingly, 
the results obtained with this equation appeared to be 
qualitatively correct. 

b. Stepwise oscillating force. Here we will briefly 
discuss the model of stepwise oscillating force to check 
the results obtained in ref. [f| with the equation which 
is not quite correct, in general. The model is defined 
as z(t) — zq(— l)[ 2 */ T ol, where To is the oscillation period 
and [x] denotes the integer part of x. It can also be rep- 
resented as 

00 

z(t) = z c (t) = z ne iW \ (35) 

n— — 00 

where u> = 2tt/t and z n are given by: z 2n = 0, z 2n +i = 
-2i/[7r(2n + 1)] with z_„ = z*. 

As mentioned above, the case of z c (t) represented as 
a superposition of harmonically oscillating functions can 
be described by assuming z c (i)-modulation to result from 
dynamical motion in a harmonic potential in the multi- 
dimensional space z = (zx,v Zl ; . . . ; z„, v Zn ; . . .). In this 
case formulas (|27p - (j2"5)l predict the same expressions (|3"Tjl - 
(j32[) for the moments. Evaluation with these expressions 

yields: mi(i) * ~°° 7i(a)(2//a; Q ), and 

Mt)^ 72(a) f(t/wr. (36) 

Here the functions 71(a) and 72(a) are defined as: 

71(a) = C(l+a)(2-2- a )sin(7ra/2), (37) 
72(a) = C(2 + a)(4-2-«)/[4r(l + a)], (38) 

where Ci x ) is the Riemann's zeta-function. 

Obtained results agree with those of ref. [6] thus con- 
firming the correctness of the method proposed in this 
work as applied to the model of stepwise oscillating force. 



c. Fluctuating force. Of great interest is also the 
case of fluctuating force F(t), i.e. fluctuating z(t). This 
case is conveniently analyzed using eqs. (|2"5|) and (f2"6")) for 
the Laplace transforms. 

For brevity, we will analyze only the long time limit. 
In addition, for simplicity, we will assume that the initial 
condition <Ji(z) in {z}-space is equilibrium i.e. Lai = 0, 
and mi(0) = (z) CTi — 0. In such a case the first moment 
m\(t) = mi(0) = so that the only value to be analyzed 
is /U2(e)' At small e — > this term can be estimated as 

Me) « f\-^{e){e.\z^{Sl)z)\e s ), (39) 

where \e z ) and (e z \ are the equilibrium states of the op- 
erator L (in this bra-ket notation (e z \ = Jdz Q and 
Wi) = IO) therefore 

Ma (t) t ~ 00 r- 1 (l + a)/ 2 (ft) a (40) 

with / = D a FQyJ (z 2 ) and the parameter f is expressed 
in terms of the correlation function K(t) — (zz(t)) = 
(e z \ze- Lt z\e z ): 

poo 

f a = dtK(t)(f)(t)/{z 2 ). (41) 
^0 

It is seen that the characteristic features of /J, 2 (t) are 
similar to those obtained in the case of oscillating force 
with / and u> replaced by / and f , respectively. Note- 
worthy is, however, that unlike this case, for fluctuating 
force n 2 (t)/t a is finite as a — » 1. 

V. CONCLUDING REMARKS. 

We have analyzed the response of CTRWs on time- 
dependent field using the rigorous method based on the 
Markovian representations of CTRW and the modulated 
field. This method is applied to describing the field effect 
on subdiffusive motion. 

Obtained formulas ((33|) -(|40 |) clearly demonstrate some 
specific features of the response of anomalous subdiffusive 
systems: 1) in the case of oscillatory force modulation the 
first moment (average displacement) is, in general, non- 
zero (even in the long time limit) and depends on the 
oscillation phase, 2) the modulated force results in the 
anomalously strong contribution to the second moment 
(dispersion) growing in time. 

In conclusion, it is worth noting that the proposed 
Markovian SLE approach ((9]) for describing the influ- 
ence of modulated external fields is applied to only one 
particular problem of the theory of force induced ef- 
fects in stochastic systems. This approach is, however, 
fairly general and can be very suitable in studying many 
time-dependent-field affected stochastic processes |l|, H, i| 
since it reduces the study to the analysis of characteris- 
tic features of time-independent operators (their spectra, 
eigenfunctions, etc.). Some of applications of the SLE 
approach ^ are currently under consideration. 
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